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ABSTRACT 

We introduce our project on galaxy evolution in the environment of rich clusters aiming at disentangling the 
importance of specific interaction and galaxy transformation processes from the hierarchical evolution of galaxies 
in the field. Emphasis is laid on the examination of the internal kinematics of disk galaxies through spatially 
resolved MOS spectroscopy with FORS at the VLT. First results are presented for the clusters MS 1008.1-1224 
(z = 0.30), CI 0303+1706 (z = 0.42), and CI 0413-6559 (F1557. 19TC) (z = 0.51). Out of 30 cluster members with 
emission-lines, 13 galaxies exhibit a rotation curve of the universal form rising in the inner region and passing 
over into a flat part. The other members have either intrinsically peculiar kinematics (4), or too strong geometric 
distortions (9) or too low S/N (4 galaxies) for a reliable classification of their velocity profiles. The 13 cluster 
galaxies for which a maximum rotation velocity could be derived are distributed in the Tully-Fisher diagram 
very similar to field galaxies from the FORS Deep Field that have corresponding redshifts and do not show any 
significant luminosity evolution with respect to local samples. The same is true for seven galaxies observed in the 
cluster fields that turned out not to be members. The mass-to-light ratios of the 13 TF cluster spirals cover the 
same range as the distant field population indicating that their stellar populations were not dramatically changed 
by possible clusterspecific interaction phenomena. The cluster members with distorted kinematics may be subject 
to interaction processes but it is impossible to determine whether these processes also lead to changes in the overall 
luminosity of their stellar populations. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: spiral — galaxies: 
clusters: individual: (MS 1008. 1-1224, CI 0303+1706, CI 0413-6559) 



1. INTRODUCTION 

Galaxy clusters provide a special environment for their mem- 
bers. In contrast to the field, the number volume density of 
galaxies is high and the relative velocities are large. The gravi- 
tational potential of a cluster is filled by the intracluster medium 
(ICM), a hot X-ray emitting gas, and the overall mass-to-light 
ratio is much larger than for the individual galaxies indicating 
the presence of vast amounts of dark matter This environment 
exerts a strong influence on the evolution of the cluster galaxies 
superposed on the (field) evolution that arises from the hierar- 
chical growth of objects and the declining starformation rates 
over cosmic epochs. Besides tidal interactions between galax- 
ies including merging as can also be observed in the field, clus- 
ter members are affected by clusterspecific phenomena related 
to the ICM (like ram-pressure stripping) or the structure of the 
cluster (like harassment). For a recent overview see third vol- 
ume of the Carnegie Observatories Astrophysics Series^. Im- 
prints of these interactions can not only be seen in present-day 
clusters, but also manifest themselves in a strong evolution of 
the population of cluster galaxies. One example is the photo- 
metric Butcher-Oemler effect of an increasing fraction of blue 
galaxies with redshift (e.g. Butcher and Oemler Jr 1978) imply- 
ing a rising percentage of starforming galaxies. Another exam- 
ple is the rapid decline of the abundance of lenticular galaxies 



(SO) from the dominant population in local clusters to a few per- 
cent at a lookback time of '-^ 5 Gyrs (e.g. Dressier et al. 1997). 
These observations have led to the question whether field spi- 
rals falling into a cluster can be subject to such morphological 
transformations that they appear as SO galaxies today. 

Independent from whether this overall scenario is true or not, 
the observed tidal interactions (either between galaxies or with 
the cluster potential) may cause substantial distortions both on 
the structure and the kinematics of the galaxies involved. In- 
deed, Rubin et al. (1999), for example, found that half of their 
sample of 89 disk galaxies in the Virgo cluster exhibit kinematic 
disturbances ranging from modest (e.g. asymmetric) to severe 
(e.g. truncated curves) pecuharities. On the other hand, many 
local cluster galaxies, for which only Hi velocity widths (in 
contrast to spatially resolved velocity profiles) were measured, 
follow a tight Tully-Fisher relation similar to field spirals (e.g. 
Giovanelli et al. 1997). This TFR connects the luminosity of the 
stellar population of a galaxy to its internal kinematics which 
are dominated by the presence of a dark matter halo (TuUy and 
Fisher 1977). 

But it is not yet clear whether the halo of dark matter and, 
therefore, the total mass of a galaxy can also be affected by 
certain interaction phenomena. In numerical simulations of the 
evolution of substructure in clusters by Spring el et al. (2001), 
the dark matter halo of a galaxy that falls into the cluster is 



' Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile (ESO Nos. 64.O-0158, 64.O-0152 & 66.A-0547) 
^ http://www.ociw.edu/ociw/symposia/series/symposium3/proceedings.html 
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truncated via tidal interactions so that the mass-to-light ratio 
of a galaxy gets reduced during its passage to the cluster core. 
Gnedin (2003) simulates the tidal field along galactic orbits in 
hierarchically growing clusters, and finds that about 40% of the 
dark halo of a massive galaxy (Vmax = 250km/s) is lost between 
z = 5 and z = 0. But the rotational velocity at ~ 5 disk scale 
lengths is predicted to hardly change (decrease by ^ 2%). 

2. OUR PROJECT 

Since in models of hierarchically growing structure clusters 
are still in the process of forming at z < 1 in the concordance 
cosmology, a higher infall rate and more interactions are ex- 
pected at redshifts 0.3 < z < 1 (e.g. Kodama and Bower 2001). 
With the availability of large telescopes, it is now feasible to 
conduct spatially resolved spectroscopy of the faint galaxies at 
these redshifts to observationally test these predictions. There- 
fore, we have performed a large campaign at the VLT target- 
ing seven distant rich clusters with 0.3 < z < 0.6. The clusters 
were chosen from a very limited list with existing HSTAVFPC2 
imaging (mainly the core regions) at the time the project started 
(1999) and that are accessible with the VLT. MS 1008-12 has 
no HST imaging but was included since it was imaged exten- 
sively during FORS science verification time. The main goal is 
to derive the two-dimensional internal kinematics of disk galax- 
ies from emission lines. In combination with measurements 
of starformation rates, luminosities and structural parameters 
we aim at disentangUng the effects of different interaction pro- 
cesses and find out about their respective actual effectiveness 
and importance for galaxy evolution. 

In this Letter, we present the results for the first three clusters 
of our survey: MS 1008-12 (z = 0.30), C10303-H7 (z = 0.42), 
and CI 0413-65 (z = 0.51). They were observed with FORSl 
in MOS mode, while the other four clusters (CI 0016+1609, 
MS 0451. 6-0305, ZwCl 1447.2+2619 & MS 2137.3-2353) had 
MXU spectroscopic observations with FORS2 requiring differ- 
ent reduction techniques that will be presented in future papers. 
While in Paper II (Jager et al. 2003) we give all the data that can 
be deduced from each single spectrum, we analyze in this Let- 
ter only the late-type galaxies which exhibit spatially resolved 
emission. One setup in the MOS mode of FORSl provides 19 
individual slitlets. For each cluster two setups have been de- 
signed with different rotation angles of the instrument. Using 
grism 600R and slitwidths of 1", the spectra have a dispersion 
of ~ 1.08 A/pixel, a spectral resolution of /? » 1200 and typi- 
cal wavelengths of A A ss 5200-7400A. In standard configura- 
tion FORS has a field of view of 6.8' x 6.8' with a spatial scale 
of 0.2"/pixel. To achieve our signal-to-noise requirements of 
S/N > 5 in the emission Unes of an /? < 23 galaxy, the total 
integration time was set to ~2hrs. Seeing conditions ranged 
between 0.7 and 1.3" FWHM. 

Spatially resolved velocity profiles Vrot(r) were determined 
from either the [O ll]3727, H/3 or [O lll]5007 emission fine. In 
eight cases, two Unes with sufficient S/N were visible which 
then were treated separately yielding consistent results. The 
spectral profile of an emission line was fitted by a Gaussian 
(in case of [O ii] by a double Gaussian) after applying a me- 
dian filter window of typically 0.6" to enhance the step- 
ping along the spatial axis. Since the apparent disk sizes of 
spirals at intermediate redshifts are only slightly larger than the 
slitwidth (1"), the slit covers a substantial fraction of the two- 
dimensional velocity field. Therefore, the spectroscopy is an 
integration perpendicular to the slit's spatial axis. Because of 



this effect, the maximum rotation velocity Vmax cannot be de- 
termined "straightforward" from the observed rotation. As de- 
scribed in Ziegler et al. (2002) and Bohm et al. (2003), we over- 
come this problem by simulating such longslit spectroscopy of 
each galaxy individually. In short, a two-dimensional velocity 
field is created assuming a specific rotation law that is weighted 
by the galaxy's luminosity profile and convolved to match the 
seeing at the time of our observations. Taking into account the 
galaxy's inclination, position angle, and disk scalelength, a syn- 
thetic rotation curve is generated with Vmax as the only remain- 
ing free parameter. Vmax is then determined by matching the 
synthetic to the observed velocity profile. For the galaxies an- 
alyzed here, we used for the model rotation curve a simple pa- 
rameterization with a linearily rising inner part that turns over 
into a flat outer part. But as is demonstrated by Bohm et al. 
(2003), Vmax is hardly changed when the Universal Rotation 
Curve by Persic et al. (1996) is used instead. The inclinations, 
position angles and scalelengths were derived via 2-D x^-fits 
of exponential disks to the galaxies' profiles in FORS images. 
The FWHM was 0.59" (MS 1008-12), 0.77" (CI 0413-65) and 
1 .0" (CI 0303+17), respectively. The fits accounted for the PSF. 
Spirals with too low an inclination (; < 20°) have not been used 
for the Vmax derivation. 

Galaxy luminosities were derived from total magnitudes of 
FORS images in the V (MS 1008-12) or / (CI 0413-65 and 
CI 0303+17) band, respectively, as measured with SExtractor 
(Bertin and Arnouts 1996). Observed magnitudes were cor- 
rected for Galactic (from Schlegel et al. (1998)) and intrin- 
sic extinction (following Tully and Fouque (1985)) and trans- 
formed to restframe Johnson B according to their spectropho- 
tometric type using model SEDs corresponding to Sa, Sb, Sc 
and Sdm, and calculated for a flat = 0.7 cosmology (Hq = 
70 km s"' Mpc"'). The average overall error in the photometry 
is estimated to be < 0.2*". 

3. KINEMATICS OF CLUSTER SPIRALS 

The MOS mode of FORSl, which was the only multiplex 
technique available for the first observations of our campaign, 
has some disadvantages for the spectroscopy of spiral galaxies 
in clusters, which are removed by the MXU mode of FORS2 
that we could use for the other four clusters of our survey. 
Firstly, the 19 slitlets of fixed length (^ 22") have only one de- 
gree of freedom for placing a shtlet onto an object. This leads 
to somewhat inefficient coverage of cluster member candidates 
and many slitlets must be filled by strongly relaxing the ideal 
selection criteria. Secondly, once a certain rotation of the field 
is chosen, all the 19 slitlets have the same orientation on the 
sky. Ideally, the sUt should be placed along the major axis of 
a galaxy to probe its rotation around the center. Although we 
have targeted each cluster with two different setups, the devi- 
ation S between slit angle and position angle was rather large 
in some cases leading to geometric distortions of the observed 
velocity profile that could not be corrected for. These galaxies 
are not used in the further analysis of the internal kinematics 
here but are valuable ingredients for future studies of e.g. star- 
formation rates or structural properties of cluster members. 

As is specified in more detail in Paper 11, our method of se- 
lecting cluster spiral candidates was different for each clus- 
ter due to the limited information of published studies. The 
most comprehensive source was available for MS 1008-12, for 
which we exploited a catalog of ~ 80 cluster members with 
published spectral types (Yee et al. 1998). A list of candidates 
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Fig. 1 . — Spatially resolved velocity profiles for 12 disk galaxies in the cluster MS 1008-12 ■dlz = 0.30. Note, panels c&d, e&f, k&l and m&n represent the same galaxy, but observed 
with different slit positions (the rotation angles of the two MOS masks differed by 67°). Panels a-p are ordered according to the distance of the galaxies from the cluster center. Values 
of the maximum rotation velocity could be determined for cases a, d, e, f, g, i, o & p indicated by the label 'TF". These galaxies enter the TFR shown in Fig. 2. In all other cases, the 
kinematics are too disturbed (for c & k this may be due to a large mismatch angle S). For the 'TF" cases, the solid lines denote the projected model curves that were produced from the 
best fitting 2-D velocity fields simulating the observational circumstances. Note, that the maximum values of these fits are not the intrinsic Vmax values needed for Fig. 2. 



in CI 0413-65 was prepared by comparing optical-nearinfrared 
colors by Stanford et al. (2002)^ to evolutionary stellar pop- 
ulation models of an updated version of Bruzual and Chariot 
(1993). For CI 0303+17, we mainly utilized a spectroscopic 
catalog of Dressier and Gunn (1992). MOS slitlets that could 
not be filled by a galaxy from our input lists were placed on 
objects selected according to their structural appearance and 
magnitudes as measured on our FORS images. If no suitable 
candidate for a late-type galaxy was available, the slitlet was 
placed on an elliptical candidate since understanding the phe- 
nomenon of galaxy transformation requires the analysis of the 
whole galaxy population of a cluster. 

Redshifts and spectral types could be determined for 12/13 
spirayelliptical members of MS 1008-12 (plus 4/2 S/E field 
galaxies), 8/1 spiral/elliptical members of CI 0413-65 (plus 
11/6 S/E field galaxies), and 10/7 spiral/elliptical members of 
CI 0303+17 (plus 15/0 S/E field galaxies). While the early-type 
galaxies will be discussed in a future paper, we here concentrate 
on the emission-line galaxies. As was pointed out by Verheijen 
(2001) only galaxies with a rotation curve that rises in the inner 
region and then clearly turns into a flat part should be used for 
a TuUy-Fisher diagram. In such a case, the measured Vmax is 
representative for the influence of the dark matter halo on the 
galaxy's kinematics and is indicative for the total dynamical 
mass of the galaxy. For our three clusters, we were able to de- 
termine Vmax for 7/5/1 different member galaxies, respectively, 
(and 1/5/1 field spirals). The remaining cluster spirals either 
have too low S/N of their emission lines to spatially analyse 



the internal kinematics (in 0/0/4 cases), have too low incUna- 
tion i or too large mismatch angle 6 (in 1/3/5 cases), or exhibit 
intrinsic distortions (in 4/0/0 cases). 

In Fig. 1 we show as an example position- velocity diagrams 
for all member galaxies in MS 1008-12 that have sufficient 
S/N. Of the 12 members, of which four were observed twice 
with different slit angles with respect to their major axis, seven 
exhibit the "classical" rotation curve shape rising in the inner 
part and turning over to a flat regime (labeled "TF"). Four mem- 
bers clearly show disturbed kinematics (panels b, j, 1, & m). The 
distortion seen in panel m (& n) most probably arises from a bar 
which is readily visible in the direct image. Two "double-hits" 
were observed with very big mismatch angles (panels c & k) so 
that their velocity profiles look peculiar. In one case (panel h) 
the object fell onto the slitlet by chance and is too weak (small) 
to derive its structural parameters. 

In the figure, the galaxies are ordered according to their pro- 
jected distance to the cluster center (ranging from 285 kpc to 
1230 kpc). There is no trend visible of the RC form as a func- 
tion of clustercentric distance, i. e. distorted velocities are 
not uniquely tied to the central region. But since our observa- 
tions cover only the region within the virial radius of the cluster 
(~ 1880 kpc), this is in accordance with dynamical models in 
which the galaxy population of a cluster is well-mixed within 
that region (Balogh et al. 2000). In particular, we most probably 
do not have any new arrivals from the field in our sample. 

4. THE TULLY-FISHER RELATION 



^ kindly provided to us by A. Stanford well in advance of their actual publication 
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In Fig. 2, we present the Tully-Fisher diagram for the distant 
cluster galaxies. Only those galaxies for which Vmax could be 
determined enter the plot. The ordinate gives the restframe ab- 
solute Johnson B-band magnitudes. We also show the position 
of those field galaxies that were serendipitously observed in the 
cluster fields and that have 0.4 <z < 0.6 (3 spirals at z = 0.61 
may be members of a background cluster behind CI 0413-65). 
For comparison, we include the distant field galaxies from the 
FORS Deep Field (Heidt et al. 2003) which were observed with 
exactly the same instrumental setup (Ziegler et al. 2002; Bohm 
et al. 2003). The linear bisector fit to this sample is flatter than 
the slope fitted to the local sample of Pierce and TuUy (1992, 
PT92), which is given with its ±3(t deviations. 

The distant cluster spirals are distributed very much alike the 
field population that covers similar cosmic epochs (0.1 < z < 1 
with the bulk of galaxies at z « 0.3 and « 0.6). No significant 
deviation from the distant field TFR is visible and the cluster 
sample has not any increased scatter, but the low number of 
cluster members prohibits any quantitative statistical analysis. 
Nevertheless, we can conclude that the mass-to-Ught ratios of 
the observed distant cluster spirals cover the same range as the 
distant field population indicating that their stellar populations 
were not dramatically changed by possible clusterspecific in- 
teraction phenomena. In particular, we do not detect any sig- 
nificant overluminosities as would be expected in the B band 
if strong starbursts had occured in the recent past of the exam- 
ined cluster galaxies. With respect to the Tully-Fisher relation 
obeyed by local galaxies (e.g. PT92), our cluster sample fol- 
lows the same trend as the FDF galaxies. Since mostly only the 
bright galaxies made it into the TF diagram, the cluster mem- 
bers occupy a region where no significant luminosity evolution 
is visible. 
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Fig. 2. — Tully-Fisher diagram of cluster spirals in MS 1008-12 (filled triangles), in 
C104I3-65 (filled squares), and C10303+17 (filled diamond). Also shown are seven 
field objects (open symbols) that were serendipitously observed. In comparison to the 
FoRS Deep Field sample of 77 field galaxies (Ziegler et al. 2002; Bohm et al. 2003) with 
0.1 < z < 1 (small open circles), the cluster galaxies are similarly distributed and do not 
deviate significantly from the linear fit to the FDF sample (solid line). The cluster mem- 
bers follow the same trend with respect to the local TFR (the fit ±3(T to the Pierce and 
TuUy (1992) field sample is given) as the distant field galaxies with the brightest galaxies 
exhibiting the smallest limiinosity evolution. 



But we emphasize that this conclusion is true only for those 
objects that enter the TF diagram. Since a significant fraction 
of our cluster galaxies can not be used for a TF analysis due 
to their distorted kinematics the above conclusions are not gen- 
erally valid for the whole cluster sample. The objects with in- 
trinsically pecuhar velocity curves may actually be subject to 
ongoing or may have experienced recent interactions. Such pro- 
cesses most probably also influence the stellar populations of a 
galaxy changing its integrated luminosity as well. Tidal interac- 
tions for example could distort the spiral arms while inducing 
starbursts at the same time leading to enhanced luminosities. 
Since we do not know where the galaxies with peculiar kine- 
matics lie in the TF plane, it is not possible to decide whether 
a particular galaxy has an increased or decreased luminosity. 
Overall, these galaxies span a very similar range in apparent 
magnitudes to the TF cluster members. 

Exploring the spatially resolved kinematics of disk galaxies 
in distant clusters has become feasible only recently. Milvang- 
Jensen et al. (2003, MJ03) studied a sample of 7 spiral galaxies 
in the cluster MS 1054.4-0321 at z = 0.83 with FORS at the 
VLT in a similar configuration as our own spectroscopy. Com- 
pared to a number of field spirals at corresponding redshifts, 
which were observed at the same time, they find that the cluster 
members have brighter B luminosities by ~ 0.5 - 1"" (~ 1 .5 - 2(7 
significance) for their rotational velocities. The difference to 
the average brightening of our cluster members is hardly signif- 
icant and may be due to a combination of low-number statistics 
and systematic deviations. But we also can not rule out that the 
differences are real and may be connected to the higher redshift 
of MS 1054 or other characteristics of that cluster. 

Metevier (2003) et al. examined galaxies in the cluster 
CI 0024+1654 (z = 0.40) with the Keck lOm-telescope. The ten 
galaxies that appear in their TF diagram have a larger scatter 
than the local PT92 sample, but show no evidence for an evolu- 
tion of the zero-point. The authors argue that processes acting 
on the cluster galaxies some time before the lookback time of 
the observations involving either starbursts or a truncation of 
star formation may have caused a decreased or increased mass- 
to-light ratio, respectively, which is manifested in the increased 
scatter. 

In a future paper we will present our analysis for the other 
four clusters of our campaign. With more cluster member 
galaxies both in the TF diagram and those with peculiar kine- 
matics, we may hopefully be able to quantitatively investigate 
the galaxy evolution in rich clusters and to give also statistical 
tests. 
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